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. Abstract. We present a study of the integrated physical properties of a sample of 28 X-ray galaxy clusters observed with 

£ — ' Chandra at a redshift between 0.4 and 1.3. In particular, we have twelve objects in the redshift range 0.4—0.6, five between 

0.6 and 0.8, seven between 0.8 and 1 and four at z >1.0, compounding the largest sample available for such a study. We 
focus particularly on the properties and evolution of the X-ray scaling laws. We fit both a single and a double 0— model with 
the former which provides a good representation of the observed surface brightness profiles, indicating that these clusters do 
not show any significant excess in their central brightness. By using the best-fit parameters of the /3— model together with the 
measured emission-weighted temperature (in the range 3-1 1 keV), we recover gas luminosity, gas mass and total gravitating 
mass out to Rsoo- We observe scaling relations steeper than expected from the self-similar model by a significant (> 3cr) 
amount in the L — T and M gas — T relations and by a marginal value in the Mtot — T and L — Mtot relations. The degree of 
evolution of the Mtot — T relation is found to be consistent with the expectation based on the hydrostatic equilibrium for gas 
within virialized dark matter halos. We detect hints of negative evolution in the L — T, A/ gas — T and L — M to t relations, thus 
suggesting that systems at higher redshift have lower X-ray luminosity and gas mass for fixed temperature. In particular, when 
the 16 clusters at z > 0.6 are considered, the evolution becomes more evident and its power-law modelization is a statistically 
good description of the data. In this subsample, we also find significant evidence for positive evolution, such as (1 + z) ' 3 , in 
the E^ 3 S — T relation, where the entropy S is defined as T/n%s and is measured at 0.1 ifeoo- Such results point toward a 
scenario in which a relatively lower gas density is present in high-redshift objects, thus implying a suppressed X-ray emission, 
a smaller amount of gas mass and a higher entropy level. This represents a non-trivial constraint for models aiming at explaining 
the thermal history of the intra-cluster medium out to the highest redshift reached so far. 



Key words, galaxies: cluster: general - galaxies: fundamental parameters - intergalactic medium - X-ray: galaxies - cosmol- 
ogy: observations - dark matter. 



1. Introduction 

The physics of the intracluster medium (ICM) is mainly driven 
by the infall of the cosmic baryons trapped in the deep grav- 
itational potential of the cluster dark matter halo. Through a 
hierarchical formation that starts from the primordial density 
fluctuations and generates the largest virialized structures via 
gravitational collapse and merging, galaxy clusters maintain 
similar properties when they are rescaled with respect to gravi- 
tational mass and epoch of formation. The shock-heated X-ray 
emitting ICM accounts for most of the baryons collapsed in the 
cluster potential (e.g. White et al. 1993). Its physical properties, 
like density and temperature, relate in a predictable way in this 
simple self-similar scenario (e.g. Kaiser 1986, 1991, Evrard 
& Henry 1991). Under the assumptions that the smooth and 



spherically symmetric ICM is heated only by gravitational pro- 
cesses (adiabatic compression during DM collapse and shock 
heating from supersonic accretion), obeys hydrostatic equilib- 
rium within the underlying dark matter potential and emits 
mainly by bremsstrahlung, one derives the following scaling 
relations between the observed (bolometric luminosity Lboi, 
and temperature T gas ) and derived (gas entropy S, gas mass 
M gas , and total gravitating mass M tot ) quantities: 

- Et /3 S oc A~ 2/3 T gas 

- E- 1 L hol ex 2 T g 2 as 

- E z M tot oc A~ 1/2 T g 3 / S 2 

- E- 1 L bol oc A 7 J 6 (£ z M tot ) 4 / 3 

- E z A:/g as oc A z -^"gas ■ 
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The factors that indicate the dependence on the evolution 
of the Hubble constant at redshift z, E z — H z /H = 
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Fig. 1. Correlations between the best-fit parameters of the 
/?— model: higher values of the core radius, r c , correspond to 
higher estimates of (3. No redshift-segregation is evident. 

1 /2 

[f2 m (l + z) 3 + 1 — n m ] ~ (for a flat cosmology with mat- 
ter density O ra ), and on the overdensity A z account for the 
fact that all the quantities are estimated at a given overdensity 
Az with respect to the critical density estimated at redshift z, 
Pc , z = 3iT| /(87rG). 

Hydrodynamical simulations (e.g. Evrard, Metzler, 
Navarro 1996, Bryan & Norman 1998, Thomas et al. 2001, 
Bialek et al. 2001, Borgani et al. 2002 and references therein) 
and observational analyses (e.g. from Mushotzky 1984, Edge 
& Stewart 1991 to the more recent work of Allen & Fabian 
1998, Markevitch 1998, Arnaud & Evrard 1999, Nevalainen 
et al. 2000, Finoguenov et al. 2001, Ettori et al. 2002) of 
the best-studied correlations between X-ray luminosity, total 
gravitating mass and gas temperature show that significant 
deviations exist between observations and the expectations 
based on self-similar scaling: (i) a difference of the order of 
30—40 per cent in the normalization of the M to t ~ T gas relation 
(Horner et al. 1999, Nevalainen et al. 2000); (ii) steeper slopes 
of the Lhoi-T, M gas -T and, possibly, Mtot-T relations, in 
particular when low temperature (Tgas < 3 keV) systems are 
also considered; (iii) excess entropy in the central regions of 
clusters and groups, with respect to the S oc T expectation 
(e.g., Ponman et al. 1999, 2003). Such deviations are currently 
interpreted as evidence for non-gravitational processes, such 
as extra heating and radiative cooling, that affect the assembly 
of the baryons in the cluster potential well (e.g. Cavaliere, 
Menci & Tozzi 1999, Bower et al. 2001, Tozzi & Norman 
2001, Babul et al. 2002, Voit et al. 2002). 

The same processes that determine the shape of the local 
relations should have affected also the evolution of the scal- 
ing laws with the cosmic epoch. However, due to the objective 



difficulty to assemble a large dataset of high redshift objects, 
there has been very little work done until now, and with contra- 
dictory results, on the observed evolution of the scaling laws. 
For example, for different subsamples of our dataset, Holden 
et al. (2002) observe no evolution of the L\, \ — T gas relation 
in a ACDM universe, while Vikhlinin et al. (2002) claim to 
have evidence at 8er level for a positive evolution, i.e. systems 
at higher redshift have higher luminosity at fixed temperature. 
In principle, any observed evolution in the scaling laws allows 
to constrain different scenarios of heating schemes and cooling 
efficiency. 

Furthermore, understanding the evolution of scaling re- 
lations involving luminosity, temperature and total collapsed 
mass, is of vital importance if the evolution of the cluster pop- 
ulation has to be used as a proxy to the cosmic evolution and, 
therefore, for the determination of cosmological parameters 
(e.g., the pioneering work of Henry & Arnaud 199 1 ; see Rosati, 
Borgani & Norman 2002, for a recent review and references 
therein). In fact, while the mass function of galaxy clusters and 
its evolution are sensitive probes of the underlying cosmologi- 
cal scenario, mass is usually inferred from observational quan- 
tities, such as the X-ray luminosity and the gas temperature, 
which are affected by non-gravitational processes. 

In this paper, we focus on the general behaviour of the 
scaling laws for X-ray clusters at high redshift, with the main 
aim of quantifying their evolution with respect to what is ob- 
served for nearby clusters. To do this, we have assembled a 
set of 28 clusters in the redshift range 0.4-1.3 observed at arc- 
sec angular resolution with Chandra (Weisskopf et al. 2000). 
This allows us to resolve on scales of few tens of kpc the 
surface brightness (i.e. the gas density) of these systems, and 
to determine a single emission weighted temperature. In the 
next section, we describe our dataset and how the physical 
quantities are measured from the observations. We present in 
Section 3 the results on the scaling laws under examination 
and in Section 4 the constraints on the observed evolution of 
these relations. We summarize our conclusions in Section 5. 
All quoted errors are at ler level (68.3 per cent level of confi- 
dence and A\ 2 — 1 for one interesting parameter) unless oth- 
erwise stated. The cosmological parameters Ho = 70hj^ km 
s _1 Mpc -1 and O m = 1 — = 0.3 are assumed hereafter as 
suggested from WMAP constraints on the anisotropies in the 
cosmic microwave background (Spergel et al. 2003). 

2. The dataset 

As of Summer 2003, we consider all the clusters available 
in the Chandra archive with redshift larger than 0.4. These 
28 galaxy clusters have redshift between 0.4 and 1.3 (median 
z — 0.73), emission-weighted temperature in the range 3-1 1 
keV (median value of 6.0 keV) and luminosity between 10 44 
erg s _1 and 1.2 x 10 46 erg s~ 4 (median value: 6.9 x 10 44 
erg s _1 ). In particular, twelve of these objects were selected 
from the ROSAT Deep Cluster Survey (RDCS; Rosati et al. 
1998, Stanford et al. 2001, Holden et al. 2002, Rosati et al. 
2003) of which we present the complete Chandra follow-up 
of objects with z > 0.8, seven from the Einstein Extended 
Medium Sensitivity Survey (MS; Gioia et al. 1990), four from 
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Table 1. Best-fit results of the spectral and spatial analysis of the sample of galaxy clusters at redshift > 0.4 and derived quantities 
at an overdensity corresponding to A z — 500 in a Einstein-de Sitter universe. A Hubble constant of 70 km s _1 Mpc -1 in a flat 
universe with fi m equals to 0.3 is assumed. Note: (i) RXJ1347-1145 has a clear hot enhancement to the South-East of the 
main X-ray emission. A sector with position angle 90° — 180° and centered on the X-ray peak has been masked in accordance 
with the analysis presented in Allen, Schmidt & Fabian (2002). (ii) In 3C295, the central active galactic nucleus and regions 
of enhanced emission associated to the radio lobes have been masked (see Allen et al. 2001a). (iii) MS1054-0321 presents 
significant substructure (e.g. Jeltema et al. 2001). The temperature and the best-fit of the surface brightness profile are estimated 
from the main body of the cluster once a circular region centered at (RA, Dec; 2000)=(10 h 56 m 55"7, — 3°37'37") and with radius 
of 36 arcsec is masked. 
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the 160 Square Degrees ROSAT Survey (Vikhlinin et al. 1998, 
Mullis et al. 2003), RXJ1113 and WGA1226 from the Wide 
Angle ROSAT Pointed Survey (WARPS; Perlman et al. 2002, 
Maughan et al. 2003; see also Cagnoni et al. 2001), RXJ1716 
is part of the North Ecliptic Pole survey (NEP; Gioia et al. 
1999, Henry et al. 2001), RXJ 1347 has been discovered in the 
ROSAT All Sky Survey (Schindler et al. 1995) and 3C295 has 
been already mapped with Einstein (Henry & Henriksen 1986). 
Subsamples of the present one have been used to constrain the 
cosmological parameters through the cluster baryonic content 
(Ettori et al. 2003) and the Iron abundance in the ICM at high 
redshift (Tozzi et al. 2003). 



The data reduction is described in details in Tozzi et al. 
(2003). We summarize here its most relevant aspects. We re- 
process the level=l events files retrieved from the archive. 
When the observations were in VFAINT mode, the pulse height 
in a 5 x 5 event island has been used. The corrections for 
Charge Transfer Inefficiency (CTI) have been applied when 
the Focal Plane temperature was 153 K. Finally, we apply 
a double correction to the ARF files. We apply the script 
apply_acisabs by Chartas and Getman to take into ac- 
count the degradation in the ACIS QE due to material ac- 
cumulated on the ACIS optical blocking filter since launch 
(see |http : / /exc . harvard . edu/| ciao/ threads/ 
applyacisabs/ and Marshall et al. 2003). This correction 
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Table 2. Mean, median and dispersion values for the best-fit 
parameters r c and (3 estimated for objects in the chosen redshift 
bins. 



range z 



h 70 kpc 



0.4 - 0.6 98, 137 (162) 

0.6 - 0.8 111,138 (76) 

0.8 - 1.0 121,166 (166) 

> 1.0 137,128 (61) 



0.615, 0.635 (0.217) 
0.639, 0.627 (0.096) 
0.542, 0.682 (0.324) 
0.808, 0.778 (0.333) 



apply to all the observation. For data taken with ACIS-I (note 
that only RXJ1347, 3C295, MS0451, MS1054 and WGA1226 
have been observed with the back-illuminated S3 CCD), we 
also apply manually a correction to the effective area consisting 
in a 7% decrement below 1 .8 keV (see Markevitch & Vikhlinin 
2001). 

2. 1. X-ray analysis 

The spectra are extracted out to a maximum radius R spe c, 
which is chosen cluster by cluster so as to optimize the 
signal-to-noise ratio in the 0.5-5 keV band where most of 
the Chandra effective area is available. They are modeled be- 
tween 0.6 and 8 keV with an absorbed optically-thin plasma 
(tbabs (mekal) in XSPEC v. 11.2, Arnaud 1996) with 
redshift fixed to the value available in literature and ab- 
sorption equal to the Galactic value (as recovered from ra- 
dio HI maps, Dickey & Lockman 1990) in correspondence 
of the X-ray peak. A local background obtained from re- 
gions of the same CCD free of any point source is used. 
The gas temperature, Tgas, metallicity and the normaliza- 
tion K of the thermal component are the only free parame- 
ters of the spectral fit performed with a Cash statistics due 
to the low number of counts per bin available (e.g. Nousek 
& Shue 1989; see also discussion on background-subtracted 
spectra at |http : / /heasarc . qsf c . nasa . qov71 docs/ 
xanadu/ xspec/ manual/ node5 7.html). 

The surface brightness profiles are obtained from the 
exposure-corrected image by fixing the number of counts per 
bin between 50 and 200, depending on the photons available, 
and is then fitted with an isothermal /?— model (Cavaliere & 
Fusco-Femiano 1976). This model provides a reasonably good 
fit to all the profiles (see plots in Appendix), apart from seven 
clusters whose measured profiles have a probability < 0.01 not 
to be a random realization of a /3-model, as drawn from a \ 2 
distribution (RXJ1416, X 2 cd = L98 ; RXJ1701 , x 2 cd _ L%; 
RXJ1525, x? od = 1-47; MS0451, X 2 cd = 1-45; MS0016, X 2 cd 
= 1.40; RXJ1113, Xrcd = 1-49; RXJ2302, x 2 cd = 1-59). This 
is not unexpected, being now well known that highly spatially 
resolved observations of nearby systems have surface bright- 
ness profiles poorly described by simple functional forms, due 
to the presence of breaks and general irregularities that cannot 
be reproduced by any smoothed curve (see, e.g., the case of 



A2390 in Allen, Ettori & Fabian 2001 and the detected front 
in RXJ1252, Rosati et al. 2003) and overall two-dimensional 
shape not strictly circular as here assumed (e.g. Maughan et al. 
2003 on RXJ1113 and Donahue et al. 2003 on MS0451). On 
the other hand, the results obtained from the application of the 
(3— model represent a robust estimates of the effective values 
under exam, such as gas and total gravitating mass, as shown 
for simulated clusters (e.g. Mohr, Mathiesen & Evrard 1999). 

We have also tested a double (3— model in the form sug- 
gested by Ettori (2000 and reference therein), which pro- 
vides a consistent physical framework for the modeling: 

S h = S ,inn(l ~ X 2 )°- 5 +^ + S (l + X 2 ) - 5 -^ = 

r c A(T)7i^ gas [/ inn i?(0.5,3Ann + l)(l-x 2 ) - 5 + 3 ^ + (1- 
f hm )B(S/3 - 0.5, 0.5)(1 + x 2 ) a5 - 3/5 ], where r c is the common 
core radius, x — r/r c , A(T) is the cooling function, no tgas is 
the central gas density and B() is the incomplete beta function. 
From a statistical point of view, none of the clusters requires to 
add two parameters (i.e. the slope /3i, m of the inner, corrected 
(3— model and the relative contribution /; nn of the central, cool 
component). The F-test (Bevington & Robinson 1992, p.205ff) 
gives probability values in the range between 0.40 (MS 1054) 
and 0.85 (MS0016), so that we can conclude that a double 
(3— model is not required by our data at the fixed threshold of 
0.95 (95 per cent confidence level). In particular, we note that 
RXJ1416 and RXJ1701, the two clusters with worst X 2 asso- 
ciated to the single j3— model, show the smaller core radii, but 
still cannot accommodate a double j3— model that provides a 
F-test statistic with probability of 0.82 and 0.67, respectively. 
On the other hand, few clusters at redshift around 0.3 not in- 
cluded in the present sample, such as MS 1008 and MS2137, 
require a double j3— model with probability of >0.99 and 0.98, 
respectively. It is worth noticing two items concerning the sur- 
face brightness of high-z clusters. First, thanks to the Chandra 
spatial resolution, we are in the condition to resolve on ~ 20 
kpc scale their inner regions. For example, we are able to de- 
fine between 4 and 7 spatial bins with a fixed number of 50 
counts within 100 kpc in the four z > 1 objects. Then, the 
evidence that their surface brightness cannot accommodate a 
double (3— model indicates that these clusters do not show any 
significant excess in their central brightness, consistent with the 
picture that systems at high redshift appear as structures in for- 
mation and not completely relaxed. We refer to a forthcoming 
paper (Ettori et al. in prep.) for a more detailed analysis of the 
distribution of the observed surface brightness, in particular in 
the central regions, as function of the object's redshift. 

We have investigated the presence of trends between the 
best-fit parameters of the j3— model, the redshift and the gas 
temperature of the objects in the overall sample by evaluating 
the significance of the Spearman's rank correlation. We obtain 
a relevant deviation (Spearman's rank correlation with a proba- 
bility P = 4 x 10~ 5 ) from the no-correlation case only for the 
well-known degeneracy between r c and (3 values, being higher 
(3 measured in correspondence of larger core radii (see Fig.^ 
for the same correlation observed in nearby and intermediate- 
redshift objects see, e.g., Ettori & Fabian 1999). A mild trend 
(P = 0.03) between r c and T gas is also observed. The latter 
correlation has a milder correspondence between j3 and T gas 
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(P = 0.17) and indicates that larger core radii are measured 
in more massive systems, with no sizable change in the outer 
slope of the density profile. As shown in Fig^and quantified 
in Table |2j the best-fit parameters r c and (3 do not show any 
significant trend with the redshift. 

The advantage of using a fj— model to parameterize the ob- 
served surface brightness is that gas density and total mass pro- 
files can be recovered analytically and expressed by simple for- 
mula: 



n gas (r) 
M tot (< r) 



«-o,gas(l + a 
3 f3 Tg &s t c 



-3/3/2 



G/lTOp 1 



(1) 



where x = r/r c , fj, is the mean molecular weight in atomic 
mass unit (=0.6), G is the gravitational constant and m p is the 
proton mass. 

In particular, the central electron density is obtained from 
the combination of the best-fit results from the spectral and 
imaging analyses as follows: 



'O.ele 



47rd ang X{l + Z) 2 X K X 10 1 

0.82 x 4?rr3 x EI 



(2) 



-14 



where d ang is the angular-size distance, K — 10" " x 
/ n p n e dV/[4:ird 2 ng (l + z) 2 ] is the normalization of the ther- 
mal spectrum measured within XSPEC and the emission in- 
tegral EI is estimated by integrating the emission from the 
spherical source out to 10 Mpc along the line of sight, EI = 
f* 1 (1 + x 2 )-^x 2 dx + Q (1 + x 2 )-^x 2 {l - cos 0)dx, with 



-Rspec/Vc and x 2 = 10Mpc/r c , 



6 = arcsin(i? spcc /r), x\ 
{fj, r c ) are the best-fit parameters of the (3— model and we as- 
sume rip = 0.82n c in the ionized intra-cluster plasma. To eval- 
uate the relative errors on the estimated quantities, we repeat 
the measurements after 1 000 random selection of a temper- 
ature, normalization K and surface brightness profile, which 
were drawn from Gaussian distributions with mean and vari- 
ance in accordance to the best-fit results. By this process, we 
obtain a distribution of the estimates of gas density, gas mass 
and total mass. The median value and the 16th and 84th per- 
centile of each distribution are adopted as central and (—1, l)er 
interval, respectively (see Table^. 

From the total mass profile, M to t(< r), obtained for each 
galaxy cluster observed at redshift z, we then evaluate Ra as 
the radius encompassing a fixed density contrast, A z , with re- 
spect to the critical density, p c z , as 



A. 



3Af t ot(< Ra) 
iirp CiZ R A 



(3) 



Observational results for nearby clusters (e.g. De Grandi 
& Molendi 2002, Pratt & Arnaud 2002) show that temper- 
ature profiles do not decline significantly (^20 per cent) 
out to an overdensity ~ 500 (for a universe with tt m —I 
and J7a =0). Therefore, to make a proper use of our single- 
temperature measurement, we estimate all the observed quan- 
tities considered in this work, i.e. total mass, gas mass and lu- 
minosity to Ra = i?5oo- Specifically, an overdensity of 500 
is assumed for an Einstein-de Sitter universe, whereas A z = 



8 


keV + 


6 


keV • 


4 


keV O 




0.0 0.2 0.4 



0.6 0. 

redshift 



1.0 1.2 



Fig. 2. Values of i?5oo as function of redshift for the cosmolog- 
ical model assumed. The conversion from temperature (mass) 
to Ra is obtained by using the best-fit results of the M — T 
relation (see Sect.3) and equation[5] As comparison of what is 
adopted in other similar work (e.g. Vikhlinin et al. 2002), we 
draw the value corresponding to a fixed radius of 2 h^Q Mpc 
(dashed line). 



500 



1 + 82 (n z - 1) / (18tt 2 ) - 39 (n z - if I (18tt 2 ) 



with n z = £l m (l + z) 3 /E z , is assumed accordingly to the 
approximation of the solution for the collapse of a spherical 
top-hat perturbation that just virialized at the observed redshift 
z (e.g. Bryan & Norman 1998). Fig.|2]shows the dependence 
of i?5oo upon the typical mass (temperature) and redshift val- 
ues spanned by the clusters in our sample. Changes of about 80 
per cent are expected for -R500 in the redshift range 0.4 and 1.3, 
whereas a fixed metric radius remains constant by definition. 

In this work, the luminosity within Ra is then estimated 
by correcting the value measured directly from the spectrum 
within i? S poc as follows 



L(< R A ) = L{< R., 



spec J 



f**(l + x 2 )- 3 f 3 x 2 dx 
EI ' 



(4) 



where EI is given in Equation |2] No further corrections on 
the measurements of the luminosity are applied, i.e. no excess 
from any centrally peaked surface brightness due to the pres- 
ence of "cooling flows" is excised. On the other hand, the lack 
of any evidence for a significant excess of the brightness in the 
core regions, as discussed above, indicates that no relevant dif- 
ferences in the total luminosity are expected when excision is 
applied (e.g. Vikhlinin et al. 2002) as we also show at the end 
of this section. 

The quoted errors on luminosity are obtained through the 
propagation of the errors obtained in XSPEC from the distri- 
bution around the best-fit values in the spectral analysis (which 
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10 

present work 

Fig. 3. Comparison between the temperatures and luminosities 
estimated in the present work and by Vikhlinin et al. (2002) 
for the 21 clusters that we have in common. {Upper panel) A 
median (mean) difference of 4 (8) per cent is present between 
our measurements of gas temperature and the values quoted 
in Vikhlinin et al. {Lower panel) The dots indicate the esti- 
mates within 2/igQ 1 Mpc, whereas the asterisks show our mea- 
surements of Lboi(< -Rsoo)- A median (mean) difference of 
(8) per cent is observed between our values and the ones in 
Vikhlinin et al. 



are the dominant contribution), and the uncertainties related to 
the j3— model in moving from i? spe c to Ra- 

In the Appendix, we show the plots of the surface bright- 
ness profiles of our clusters, along with the /3— model obtained 
through a least-squares fitting process and the relative positions 
of Ra, i? S pcc and R2 a - The latter radius indicates up to where 
the signal is above 2a with respect to the best-fit value of the 
background and gives an indication of the region of the cluster 
where a diffuse emission is detectable. 

Out of 28 clusters, 17 have i?2 CT that is smaller than -R500 
by a factor between 1.1 and 3.0 (median of 1.4) and need ex- 
trapolation to compute their quantities. The spectral aperture, 
i? spcc , corresponds roughly to an overdensity of ~ 2000 and 
is about 0.6 i?2<r and 3.6 times larger than the cluster core ra- 
dius (median values). It is worth noticing that more than 75 per 
cent of the cluster luminosity, on average, is however sampled 
within R spcc because of the rig as dependence of the X-ray ther- 
mal emission. 

Finally, we compare our estimates of T gas and Lboi with 
the values quoted by Vikhlinin et al. (2002). These authors 
have 21 clusters in common with us (including also CL0024 
that is not listed here because at redshift 0.395 and, there- 
fore, below our cut at z = 0.4). The agreement is remark- 
ably good, with a median ratio between our and Vikhlinin 



Table 3. Best-fit results on the scaling relations corrected by 
the cosmological factor E z (dashed and dotted lines in the 
plots; see equation [5}. The temperature, T gaS 6, is in unit of 
6 keV; the luminosity, Lboi,44, 
tal mass, M to t,i4, in 10 1 



in 10 13 /i 70 5 ^ 2 Mq. When the slope A is fixed, we estimate 
the error- weighted mean of (logY — AlogX) and evalu- 
ate the error after resampling Y and X by 1 000 times ac- 
cording to their uncertainties. The scatter on Y is measured 

9 I 1 / 2 

T, 3 =i.N ( l °S Yj-a-A log Xj) /N . Note that the 
scatter along the X-axis can be estimated as eri og x — 0iog y I A. 



in 10 44 /i 70 2 erg s -1 ; the to- 
^/i^M©; the gas mass, M gas ,i 3 , 



relation (Y - 


-x) 


a 


A 


ay 


E z 1 Lbol,44 




0.79(±0.07) 
1.00(±0.02) 


3.72(±0.47) 
2.00 (fixed) 


0.35 
0.41 


-Ez-Mtot.14 - 


Tgas, 6 


0.75(±0.03) 
0.74(±0.02) 


1.98(±0.30) 
1.50 (fixed) 


0.15 
0.15 


-E z M gaSi i3 - 


Tgas, 6 


0.79(±0.03) 
0.92(±0.01) 


2.37(±0.24) 
1.50 (fixed) 


0.17 
0.22 


E z 1 Lbol,44 


- E z M tot ,i4 


-0.63(±0.32) 
0.09(±0.03) 


1.88(±0.42) 
1.33 (fixed) 


0.47 
0.52 



et al. measurements of temperature of 1.04 (mean: 1.08; see 
Fig. [3} and estimates of luminosity (at 2h$Q Mpc) of 1.00 
(1.08 in average). Larger (> \2a\ when comparing our values 
vs Vikhlinin et al. ones) deviations are observed in the tem- 
perature estimates of RXJ1347 (-4,7a), RXJ1701 (-3.0ct), 
CL0024 (-2.3a), 3C295 (-2.3a) and MS 1054 (+2.5er), and 
in the luminosity measurements of 3C295 (+8.2cr), RXJ1416 
(+6.2(7-), RXJ1317 (+3.1(r) and MS0302 (+2.3a). We notice 
that all these objects (apart from MS 1054 that presents a signif- 
icant substructure that we mask out in our analysis as described 
in the caption of Table 1) are classified by Vikhlinin and collab- 
orators as systems with a "possible cooling core" (even though 
they do not describe how they quantify the sharpener of the cen- 
tral brightness). The core emission is then excised by Vikhlinin 
et al. during the process of the determination of the gas temper- 
ature and luminosity, whereas we maintain it in our procedure. 



3. The scaling relations at 0.4 < z < 1.3 

A sensible way to test the hierarchical formation of galaxy clus- 
ters and the influence of non-gravitational processes on their 
self-similarity at different epochs is to study how the scaling re- 
lations behave and evolve at high redshifts. In practice, by using 
the observed (L, T) and derived (Af to t, M gas ) physical proper- 
ties measured within a given overdensity A z with respect to 
the critical overdensity p c z , one wants to investigate whether 



the quantities E z M tot /Tgls and E z 1 L ho \/T^ s [and the com- 
plementary ones, E z M gas /Tgls and E, 1 L hol / (E z M tot ) 4 / 3 ] 
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Fig. 4. L — T relation without (left) and with {right) correction by E z . (Upper panels) Dotted line: slope fixed to the predicted 
value of 2. Dashed line: slope free. The solid lines represent the local best-fit results (from thinnest to thickest line): Markevitch 
(1998; with core excised in cooling flows clusters), Arnaud & Evrard (1999), Novicki et al. (2002; for objects at z < 0.3) and 
Novicki et al. (2002; for objects at 0.6 > z > 0.3). Solid line: Ettori et al. (2002; thinest line), Allen et al. (2001). (Bottom 
panels) Plot of the A\ 2 distribution for the one interesting parameter B(B Z ) given in eqn.0 Each solid line corresponds to a 
local scaling relation plotted with the same thickness in the upper panel and compared to all the 28 objects with z > 0.4 in our 
sample. The 2 and 3<7 limits are shown as dot-dashed lines. 



show, or not, any evolution with redshift and/or dependence 
upon one of the physical properties, typically the temperature. 

The behavior of these scaling laws are examined first in 
their normalization and slope by fitting the logarithmic relation 

\ogY = a + AlogX (5) 

between two sets of measured quantities {Xj} and {Yj}. We 
use the bisector modification (i.e. the best-fit results bisect 
those obtained from minimization in vertical and horizontal 
directions) of the linear regression algorithm in Akritas & 
Bershady (1996 and references therein, hereafter BCES) that 
takes into account both any intrinsic scatter and errors on the 
two variables considered as symmetric. The uncertainties on 
the best-fit results are obtained from 10000 bootstrap resam- 
pling. The results on the best-fit normalization and slope for the 
scaling laws here investigated are quoted in Tableland shown 



as dotted (when the slope is fixed to the predicted value) and 
dashed (when the slope is a free parameter) lines in Figures |4| 
|3|6land|7l 

In general, we observe a steeper slope than expected from 
self-similar model in all the scaling laws under exam, with 
L — T (Fig. |4ji and Af gas — T (Fig.|6jl relations being the ones 
with deviations larger than 3er (best-fit of A = 3.72 ± 0.47 
and 2.37 ± 0.24, respectively), whereas M tot - T (Fig. |3} 
and L — M tot (Fig. deviates by about l.ha from the ex- 
pected values of A = 1.5 and 1.33. This is a trend usually 
well-observed at lower redshift and that affects all the quanti- 
ties depending directly on the measurements of the amount of 
gas, like luminosity and gas mass (e.g. Edge & Stewart 1991, 
Fabian et al. 1994, Mushotzky & Scharf 1997, Allen & Fabian 
1998, Mohr et al. 1999, Novicki et al. 2001, McCarthy et al. 
2002, Ettori et al. 2002). This confirms that simple gravita- 
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tional collapse is not the only process that governs the accre- 
tion of baryons in the potential well, but some combination 
of properly efficient cooling and extra heating is required to 
take place before or during the hierarchical formation of clus- 
ters (Kaiser 1991, Evrard & Henry 1991, Cavaliere, Menci & 
Tozzi 1999, Tozzi & Norman 2001, Balogh et al. 2002, Voit et 
al. 2002, Tornatore et al. 2003). Scaling the X-ray luminosity 
/ 2 T 2 , David et al. (1993) suggested 



as L oc n 2 T 2 R 3 

gas 

that if the gas mass fraction, / g 



Mg as /Mtot, is propor- 



tional to T 5 , the steepening of the L — T relation can be ex- 
plained. However, we rule out any statistically significant de- 
pendence of / gas on the gas temperature in two ways. First, 
we measure a probability of 66 per cent that no-correlation is 
present in the distribution of the values in the plane / gas — T 
once the Spearman's rank coefficient is estimated (by fitting 
a power-law to the data, we measure a slope of 0.5 ± 1.0). 
Then, we try to constrain the best-fit parameters of the model 
L ~ T A1 (Mgas/Mtot)" 42 given the estimated gas luminosities, 
gas masses and total gravitating masses. As shown in Fig. [8] Al 
is consistent with a value of 3 and different from 2 by ~ 3er, 
whereas A2 is in agreement with the predicted value of 2 within 
la (A2 = 2A8t° Q f ). 

As originally pointed out by Horner et al. (1999) and 
Nevalainen et al. (2000), samples of nearby clusters indicate 
that systems at given temperature are systematically less mas- 
sive than the counterparts measured in numerical simulations 
(e.g. Evrard et al. 1996, Mathiesen & Evrard 2001, Borgani 
et al. 2002), where normalizations of the M to t — T relation 
are larger by about 30 per cent or more than observed, once 
the slope is fixed to the expected value of 3/2. We confirm 
this result at z > 0.4, where we measure a normalization of 
(0.37 ± 0.02) x lO^/i^M© keV- 3 / 2 , that is indeed - 30 
per cent lower than what obtained in the hydrodynamical sim- 
ulations of Mathiesen & Evrard (2001; see Fig. |5jl, but well 
in agreement with previous observational results obtained for 
cluster samples at lower redshifts (e.g. Nevalainen et al. 2000, 
Allen, Schmidt & Fabian 2001). 

We note that the observed scatter in these relations indi- 
cates a correlation between total/gas mass and gas temperature, 
which is as tight at high redshift as measured locally. This sug- 
gests that a marginal role is played by evolutionary effects, such 
as merging history and non-gravitational heating, during the as- 
sembly of baryonic and dark matter into the clusters. A larger 
scatter is present in the L— T and L— Mtot relations. This is just 
the consequence of the higher sensitivity of X-ray luminosity 
to the details of the gas distribution and how it is affected by 
both gravitational dynamics and non-gravitational processes. 



4. On the evolution of the scaling laws 

To constrain the evolution in the considered scaling laws, we fix 
(a, A) to the best-fit results obtained from a sample of objects 
observed at lower redshift, (a, A), and evaluate the confidence 
interval through a least-square minimization on the parameter 
B in the relation 



logF = a + AlogX + Slog(l + z). 



10 : 



o 



+ 0.4<z<0.6 
X0.6<z<0.8 
*0.8<z<1.0 
. z>1 .0 




1 



gas 



(keV) 
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Fig. 5. Mtot — T relation. (Upper panels) Dotted line: slope 
fixed to the predicted value of 3/2. Dashed line: slope free. 
The solid lines represent the local best-fit results (from thinnest 
to thickest line): Ettori et al. (2002), Finoguenov et al. (2001), 
Allen et al. (2001). The dash-dot line indicates the best-fit re- 
sults from hydrodynamical simulations in Mathiesen & Evrard 
(2001, Table 1, with temperatures from simulated spectral anal- 
ysis in the band 0.5-9.5 keV). (Bottom panels) Plot of the 
A\ 2 distribution for the one interesting parameter B z given 
in eqn.0 Each solid line corresponds to a local scaling relation 
plotted with the same thickness in the upper panel and com- 
pared to all the 28 objects with z > 0.4 in our sample. The 2 
and 3cr limits are shown as dot-dashed lines. 



In particular, for a given grid of values of {Bi}, we search for 
the minimum of the merit function 



E 



[log Yj-a-A log Xj - Bi log(l + Zj )] ' 

-t-2 



(7) 



(6) 



where the errors on the best-fit local values, e& and e-r, are 
considered and propagated with the uncertainties, e\ gx — 
ex/(X In 10) and e\ og Y — ey/(Y In 10), on the measured 
quantities. 
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Table 4. The best-fit values of the evolution parameter B, B z (see equations[6]and0and comments in Section 4) are obtained by 
comparing to the local best-fit results, indicated in the second column and discussed in the text, the 28, 16 and 1 1 clusters with 
redshift larger than 0.4, 0.6 and 0.8, respectively. The case of no evolution corresponds to a value of B z equal to 0. The L — T 
relation has been studied also not including the correction by E z . This relation is then expected to be proportional to E z and B 
should be measured in the range 0.6-0.9 in the adopted cosmology. 



relation 



local ref. 



> 0.4) 



> 0.6) 



> o.; 



L-T 

(fig-E) 



Markevitch 98 
Arnaud & Evrard 99 
Novicki et al. 02 (z < 0.3) 



0.62(±0.28) 
0.98(±0.20) 
0.54(±0.31) 



B 

0.04(±0.33) 
0.22(±0.25) 
-0.06(±0.37) 



0.04(±0.39) 
0.24(±0.30) 
-0.02(±0.43) 



Novicki et al. 02 (0.3 < z < 0.6) 0.10(±0.62) -0.48(±0.73) -0.34(±0.83) 









B z 




E~ X L - T 


Ettori et al. 02 


-1.04(±0.32) 


-1.48(±0.38) 


-1.46(±0.44) 


(fig.ffl 


Allen et al. 01 


-0.16(±0.23) 


-0.72(±0.28) 


-0.78(±0.32) 


EzMtot — T 


Ettori et al. 02 


-0.12(±0.39) 


-0.14(±0.46) 


-0.06(±0.52) 


(fig.0 


Finoguenov et al. 01 


-0.26(±0.17) 


-0.16(±0.21) 


-0.12(±0.26) 




Allen etal. 01 


0.02(±0.39) 


0.04(±0.45) 


0.06(±0.52) 


£ z M gas - T 


Ettori et al. 02 


-0.26(±0.23) 


-0.54(±0.27) 


-0.54(±0.32) 


(fig-U 


Mohr et al. 99 


-0.10(±0.20) 


-0.42(±0.24) 


-0.42(±0.28) 


E Z 1 L — E z Mtot 


Ettori et al. 02 


-0.60(±0.17) 


-1.10(±0.20) 


-1.12(±0.25) 


(fig-E} 


Reiprich & Bohringer 02 


0.18(±0.20) 


-0.86(±0.25) 


-0.78(±0.32) 



The local and intermediate-redshift relations considered in 
this work have been rescaled to the assumed cosmology and 
overdensity as described in the Introduction. On the other hand, 
many results on the clusters scaling laws, and particularly on 
the L — T relation, are presented in literature with the omission 
of the cosmological factor E z (that should be properly con- 
sidered only when the quantities are measured within a fixed 
overdensity with respect to the critical density as done in the 
present work). Therefore, we have considered in the follow- 
ing analysis both the inclusion and the omission of the factor 
E z to correct the quantities whose relations are investigated. In 
particular, we refer to the parameter B for the evolution in the 
scaling relations not corrected by E z and to the parameter B z 
for the evolution in the corrected scaling laws. 

It is worth noticing that the cosmological factor E z is de- 
fined as H z /H = [n m (l + z) 3 + 1 - Q m ] 1/2 (for a flat 
cosmology with matter density f2 m ) and is exactly equal to 
(1 + z) in an Einstein-de Sitter universe and proportional 
to (1 + z )°-6/o.9 in the r edshift range here considered for an 
assumed ACDM model with Q m = 0.3. 1 Furthermore, consid- 



1 The slope of the dependence of E z upon (1 + z) is the values of 
the ratio log_E z /log(l + z) for an assumed cosmology and a given 
range in redshift. 



ering its dependence upon the redshift, this factor affects sig- 
nificantly the quoted values of, e.g., mass and luminosity for 
objects at high— z. For example, the X-ray luminosity has to be 
corrected by a factor of 0.81 at z = 0.4 and 0.51 at z = 1.2, 
making these corrections fundamental to proper calibrate the 
scaling relations in a sample of clusters distributed over a wide 
interval in redshifts. Thus, the omission of the factor E z , that 
affects each quantity in a different way at different redshift as 
stated above, has to be properly considered once one wants to 
assess any evolution of the physical quantities that is not due to 
the expansion of the Universe and to the rescaling of the char- 
acteristic density. 

With this notation in hand, we should observe B z = and 
B equal to the expected dependence of the scaling law upon 
E z if no evolution of the physical quantities is present. In the 
following analysis, only the evolution of the L — T relation is 
presented by separating B and B z . One should bear in mind 
that the L — T relation is then proportional to E z implying that 
B should be measured in the range 0.6 — 0.9. 

We describe below the constraints that we obtain on the 
parameters B and B z , whose x 2 distribution is plotted in the 
bottom panel of the figure that shows the corresponding seal- 
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B z (evolution) 



Fig. 6. M gas — T relation. (Upper panels) Dotted line: slope 
fixed to the predicted value of 3/2. Dashed line: slope free. The 
solid lines represent the local best-fit results (from thinnest to 
thickest line): Ettori et al. (2002), Mohr et al. (1999). (Bottom 
panels) Plot of the A\ 2 distribution for the one interesting pa- 
rameter B z given in eqn. Each solid line corresponds to a 
local scaling relation plotted with the same thickness in the up- 
per panel and compared to all the 28 objects with z > 0.4 in 
our sample. The 2 and 3er limits are shown as dot-dashed lines. 



ing relation. All the results, obtained also by selecting different 
redshift-range in our sample, are presented in Table0] 

4. 1. Evolution in the L-T relation 

We compare the distribution of our measurements in the L — T 
plane (Fig. 0} with the best-fit results obtained for homo- 
geneous and representative sets of nearby clusters from (i) 
Markevitch (1998), who analyzed a sample of 31 clusters in 
the redshift range 0.04-0.09 with measured ASCA temperatures 
and ROSAT luminosities and excises systematically their inner 
100 kpc to avoid any contamination from central cooling cores, 
(ii) Arnaud & Evrard (1999) that include in their sample 24 
clusters with weak or no cooling flows in their cores and X-ray 




0.1 I , , 

1 1 

E z M tot (h 7 -J 10 14 M ) 




B z (evolution) 

Fig. 7. L — Mtot relation. (Upper panels) Dotted line: slope 
fixed to the predicted value of 4/3. Dashed line: slope free. 
The solid lines represent the local best-fit results (from thinnest 
to thickest line): Ettori et al. (2002), Reiprich & Bohringer 
(2002). (Bottom panels) Plot of the A\ 2 distribution for the 
one interesting parameter B z given in eqn. Each solid line 
corresponds to a local scaling relation plotted with the same 
thickness in the upper panel and compared to all the 28 objects 
with z > 0.4 in our sample. The 2 and 3er limits are shown as 
dot-dashed lines. 



measurements from GINGA, ASCA and Einstein, (iii) Novicki, 
Sornig & Henry (2002, Tab. 4 and 5), that use a sample of 53 
low (z < 0.3) and 32 intermediate (0.3 < z < 0.6) redshift 
clusters observed with ASCA . Considering that the best-fit re- 
sults obtained from these samples do not include the correction 
by the cosmological factor E z , we consider here in the same 
way our measurements without this correction (Fig. [4] panel 
on the left). In this case, the L — T relation should evolve like 
E z w (1 + z) - 6 / - 9 (fora'Un = 1-Q A = 0.3" universe) and 
we should measure B = 0.6 — 0.9, if self-similar predictions 
are correct. 

Moreover, we compare the distribution of our data cor- 
rected by E z with the cosmology-corrected relations obtained 



S. Ettori et al.: Scaling laws in high-z X-ray Galaxy Clusters 



11 



< 




o I- i i i i j 

12 3 4 5 

A1 

Fig. 8. Slopes in the relation Lbol,44 = ^40 x 
^as, 6 (%M3/Mtot,i4) i2 (l + z) A3 obtained by mini- 
mization of the x 2 (best-fit results: AO = 2.0, Al = 2.8, 
A2 = 1.9, A3 = 0.9, xldn = 29.2 with 24 d.o.f.) and 
after marginalization on the parameters AO and A3. The 
dashed lines indicate the predicted values of 2 from the 
simple gravitational collapse model assuming bremsstrahlung 
emission. 



from resolved temperature of (iv) 20 nearby T > 3 keV ob- 
jects observed with BeppoSAX (Ettori, De Grandi & Molendi 
2002) and (v) six clusters with redshift between 0.10 and 0.46 
observed with Chandra from Allen et al. (2001). A predicted 
E^ 1 L — T relation constant in time, and B z = 0, should be 
then observed. 

When we compare the distribution of our data not corrected 
by E z with the local relations obtained from samples (i)-(iii) 
(Fig. |4] panel on the left), we measure B between 0.1 and 
0.6. In particular, when we use as local reference the result 
by Arnaud & Evrard (1999), B becomes significantly positive 
(0.98 ± 0.20), but with a X 2 statistic ( X 2 = 51.6, d.o.f. = 27) 
that is worst by a Ax 2 ib 30 than the others measurements, 
mainly because of the relative small error on the slope A. These 
measurements of B are slightly lower than, but still consistent 
with, the expected value between 0.6 and 0.9. The evidence 
that B tends to be lower then the predicted value seems to 
suggest the presence of a slight negative evolution, with clus- 
ters at higher redshift having lower bolometric luminosities for 
given temperature. However, the small departure from B = 
in our L — T relation not corrected by E z agrees with previ- 
ous determinations obtained from Mushotzky & Scharf (1997), 
Fairley et al. (2000) and Holden et al. (2002). In particular, 
Holden et al. analyze a subsample of the present dataset with 
only 12 objects, nine of those were exposed with Chandra . 
On the other hand, this value of B is significantly lower than 



the result presented in Vikhlinin et al. (2002). These authors 
have assembled 22 clusters in the redshift range 0.39 — 1.26 
which were available with sufficient long exposure to ensure 
a good temperature determination in the Chandra archive as 
of Spring 2002. Single temperatures (with column density and 
metallicity fixed to Galactic value and 0.3 times the solar abun- 
dance, respectively) were estimated collecting photons within 
0.5-1 h§Q Mpc. X-ray luminosities were measured within 2 
h^Q Mpc. Both in the spectral and spatial analysis, the central 
regions with r —100 h§ kpc were excised if a sharply peaked 
surface brightness was detected. Once compared with the local 
L — T distribution in Markevitch (1998), where the same pro- 
cedure in the determination of temperatures and luminosities 
has been adopted, Vikhlinin et al. claim a detection significant 
at 8cr confidence level of a positive evolution (B = 1.5 ± 0.2). 

To check what causes this relevant difference, we re- 
evaluate B for the sample by Vikhlinin et al., using the re- 
sults by Markevitch (1998) for the low-z reference relation. 
We apply to this sample the same analysis procedure applied 
to our set of 28 clusters, which is based on minimizing the \ 2 
in equation while propagating uncertainties both in normal- 
ization and slope. We obtain B = 1.1 ± 0.3, with a deviation 
from zero of only 3.5<r. This significance decreases to 2.2<r 
(B = 0.62 ± 0.28, x 2 = 22.0 with 27 d.o.f.) when our es- 
timates of L(< -R500) (see asterisks in Fig. [3} are considered 
for our sample of 28 objects. This value of B is well in agree- 
ment with the predicted range of 0.6 — 0.9 expected when a 
face value of luminosity not corrected by E^ 1 is considered. 

Furthermore, it is worth noticing that the procedure adopted 
here defines the radius, at which all the physical properties 
of a system are recovered (see discussion in Sect. 2) as a func- 
tion of its mass (temperature), redshift and given overdensity. 
As we show in Fig. [2] the value of Ra varies considerably in 
the redshift range 0.4-1.3 with respect to a fixed metric radius. 
Therefore, systematic differences between quantities (such as 
X-ray luminosity) measured at 2 /i^ Mpc and -R500 are ex- 
pected and, in fact, are observed. If we extrapolate further our 
estimates of luminosity up to R200, on average increase of 8 
per cent in L(< -R200) with respect to L(< -R500) is measured 
and a more significant evolution (B » 0.7 ± 0.3) is observed. 

Therefore, we conclude that the mismatch between our and 
Vikhlinin et al. (2002) results is due to, in decreasing order of 
relevance, (i) the different fitting procedure (that reduces the 
significance of the evolution by propagating properly the errors 
on the best-fit parameters of the local relation), (ii) being our 
sample larger at higher redshift (e.g., we have 16 galaxy clus- 
ters at z > 0.6 and 4 at z > 1, whereas Vikhlinin et al. have 9 
and 2 objects, respectively), (iii) systematic differences in the 
definition of the reference radius and in the procedure to esti- 
mate the gas temperature and luminosity (see also discussion 
at the end of Sect. 2.1 and Fig. [3}. 

We consider now the cosmological corrected L—T relation, 
that should not evolve in the self-similar scenario. We compare 
our observed luminosities, which were corrected by the E z fac- 
tor, to the best-fit relations obtained from the samples (iv-v) 
discussed above. We measure B z = —1.04 ± 0.32 when our 
data are compared to the determination by Ettori et al. (2002) 
at median redshift of 0.05, and B z = —0.16 ± 0.23 when com- 
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paring to the results by Allen et al. (2001) for clusters at median 
redshift of 0.28 (Fig.|4l panel on the right). When we consider 
only the 16 objects with z > 0.6, B z is measured more signifi- 
cantly negative: B z = -1.48 ± 0.38 and -0.72 ± 0.28, respec- 
tively. We also note that the \ 2 is definitely improved (values 
between 2 and 9 with 15 d.o.f.) also when the B parameter is 
evaluated (B values in the range —0.5 ± 0.7 and ± 0.3). 

In agreement with a value B lower than 0.6, a B z < indi- 
cates that a marginal detection of negative evolution (i.e. clus- 
ters at higher redshift have lower luminosities for a given T) 
appears when the relation is properly corrected by its cosmo- 
logical dependence and compared to local estimates obtained 
under the same condition (i.e. all the photons are considered to 
recover both temperature and luminosity and the E z factor is 
taken into account as done in Allen et al. 2001 and Ettori et al. 
2002). 



4.2. Evolution in the M t , 



T relation 



We compare our distribution in the E z M to t — T plane with the 
local relations measured by (i) Ettori et al. (2002; BeppoSAX 
data), (ii) Finoguenov, Reiprich & Bohringer (2001) for their 
subsample of 26 clusters with spatially resolved temperatures 
and emission-weighted values larger than 3 keV as measured 
with ASCA (we have re-fitted their relation after rescaling the 
mass measurements in our cosmological model and correct- 
ing them by the corresponding factor E z ) and (iii) Allen et al. 
(2001; Chandra data). In Figure [5] we show our constraints on 
the evolution parameter B z . Overall, there is no evidence of 
any evolution with values of B z always consistent with zero: 
B z = -0.12 ± 0.39 (Ettori et al. 2002), B z = -0.26 ± 0.17 
(Finoguenov et al. 2001) and B z = 0.02 ± 0.39 (Allen et al. 
2001). Also if we restrict the measurements to the subsample 
of clusters with z > 0.6, we obtain B z consistent with zero: 
B z = -0.14 ± 0.46 (Ettori et al. 2002), B z = -0.16 ± 0.21 
(Finoguenov et al. 2001) and B z = 0.04 ± 0.45 (Allen et al. 
2001). 



4.3. Evolution in the M. 



gas 



T relation 



Fig-IUshows that there is no significant evolution in the Af gas — 
T relation corrected by E z when it is compared with local esti- 
mates from Mohr et al. (1999) and Ettori et al. (2002). For this 
relation, we find B z = -0.10 ± 0.20 and -0.26 ± 0.23 for 
these two reference analyses, thus corresponding to a detection 
of a negative evolution at < ler confidence level (in case of no 
evolution, B z — 0). A slightly more significant evolution is de- 
tected when z > 0.6 clusters are used: B z = —0.42 ± 0.24 and 
—0.54 ± 0.27, respectively. Still, it is worth noticing that, as 
observed in the L — T relation, if any hint of evolution can be 
detected, this points toward a negative evolution, i.e. clusters 
tend to have lower X-ray luminosity and gas mass at higher 
redshift for given temperature. We underline that this is only a 
marginal (< 2a) evidence, even if confirmed from both mea- 
surements of X-ray luminosity and gas mass. 
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Fig. 9. Distribution of the entropy values measured at 0.1 x 
i?200 versus gas temperature. The distribution shows a 
Spearman's rank correlation significant at 99.9 per cent against 
the null-hypothesis that no-correlation is present. A rescal- 
ing factor of E^ Z is used to compare these values to 
what is expected from hydrodynamical simulations, S rs 
50 (/ gas /0.106)- 2 / 3 T kcV /i7 4/3 kev cm 2 (solid line, from 
Ponman et al. 1999 and normalized to the median gas frac- 
tion measured in our sample; dashed line, S 
Ponman et al. 2003). 
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4.4. Evolution in the L - M tot relation 

The cosmologically corrected L — Mtot relation (Fig. [7) ap- 
pears consistent with a case of slightly negative evolution (i.e. 
B z < 0) that affects particularly the X-ray luminosities at 
redshift greater then 0.6: we measure B z — —0.60 ± 0.17 
(-1.10 ± 0.20 for objects with z > 0.6) and 0.18 ± 0.20 
(—0.86 ± 0.25 at z > 0.6) when compared with the results 
in Ettori et al. (2002) and Reiprich & Bohringer (2002), re- 
spectively. Note that we have obtained and considered the best- 
fit constraints on the E~ 1 L DO \ — E z M^q relation for the 85 
objects (out of 106) in the Reiprich & Bohringer (2002) sam- 
ple with quoted X-ray temperature larger than 3 keV. Even 
though the scatter that we measure in this relation is the highest 
among the scaling laws investigated, the trend observed for sys- 
tems with lower luminosities for given mass at higher redshifts 
agrees with the results presented on the L — T relation. 

4.5. Evolution in the S - T relation 

As a further diagnostic to trace the evolution of the ICM prop- 
erties, we analyze the relation between entropy, S, and gas tem- 
perature. Starting from the work by Ponman et al. (1999), en- 
tropy has received a considerable attention as a diagnostic of 
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the physical processes that establish the thermal status of the 
diffuse baryons in clusters (e.g., Tozzi & Norman 2001, Voit et 
al. 2002). For the first time, we attempt here to trace the evo- 
lution of the ICM entropy, whose analysis has been limited so 
far to low-z clusters. To this purpose, we measure the quantity 

5 = Tg as /n;y 3 at 0.1 x -R200 within each cluster as done by 
Ponman and collaborators for nearby systems. With the above 
definition, the quantity S relates to the thermodynamical en- 
tropy K of the system as K oc log S. 

In Figure[9] we plot the values of S as a function of the gas 
temperature. The predicted scaling of S ~ T is in agreement 
with the overall fit of the data [S = (407 ± 20) x T g ° a f ±0 20 
keV cm 2 ]. A mild negative evolution (B z = -0.14 ± 0.04, 
Xmin = 151 with 27 degrees of freedom) is observed when all 
our sample is compared with the local estimate of the "entropy 
ramp" by Ponman et al. (2003; see Fig.|9j panel on the left). 
However, the high x 2 value warns that this is not a good mod- 
elization of the data distribution, mainly for their large internal 
scatter. If we consider only the objects at z > 0.6, 0.8 and 1, 
respectively, we observe consistently B z « 0.3 with a devia- 
tion from no evolution (B z = 0) at the 2 — 3cr confidence level, 
with a reduced \ 2 of about 1 (e.g., B z — 0.26 ± 0.06, x 2 /d.o.f. 
= 18.5/15 at z > 0.6, B z = 0.32 ± 0.08, x 2 /d.o.f. = 10.6/10 
atz > 0.8). 

4.6. Discussion on the observed evolution 

As for the E~ l L — T and E z M gas — T relations, we detect 
a mild negative evolution [B z ~ (—0.1, —1.0), with a devia- 
tion from the non-evolution case, B z = 0, by less than 3.2a] 
that becomes more relevant when subsamples with higher cut 
in redshift are considered. For the 16 objects with z > 0.6, we 
measure B z in the range [—1.5, —0.4] with deviation from zero 
significant at 3.9a (see Table |4ji. This result points to a sce- 
nario in which clusters at higher redshift have in general lower 
X-ray luminosities and gas masses than what predicted from 
the self-similar model. This trend appears more clearly when 
the entropy of these hot (T gas > 3 keV) systems, evaluated at 
QAR200, is compared with local determinations. A mild posi- 
tive evolution (B z ~ 0.3) is then detected in the E^ 3 S — T 
relation, in particular when objects at redshift larger than 0.6 
are considered. 

The evolution of the L — T and M gas — T relations is a 
sensitive test for non-gravitational heating models (e.g. Tozzi 

6 Norman 2001, Bialek et al. 2001) and carries information 
about the source of this heating and the typical redshift at which 
it took place. The numerical simulations by Bialek et al. (2001) 
suggest that M gas is lower at higher redshift for fixed temper- 
ature, when the ICM is heated to an initial entropy of about 
100 keV cm 2 , well before the cluster collapse takes place. This 
effect is more evident for smaller systems, thus making the 
Mgas — T relation steeper at increasing redshift. A marginal, 
negative evolution in the L — T relation is predicted by the 
semi-analytical model by Tozzi & Norman (2001), when an 
epoch-independent entropy background of K = 0.3xl0 34 ergs 
cm 2 g -5 / 3 is combined with cooling in a ACDM universe (see 
their Fig.l 1 for emission-weighted temperature above 2 keV). 



We conclude that the observed negative evolution in the L — T 
and Mg as — T relations are in agreement with predictions from 
models that pre-heat the ICM with a non-evolving entropy floor 
at the level of about 100 keV cm 2 . In any case, this picture has 
been recently questioned by Ponman et al. (2003), who find no 
evidence for isentropic cores in groups. 

A more controversial point is whether pre-heating can also 
account for the positive evolution of entropy in central cluster 
regions. If we are looking just at the effect of an entropy floor, 
which has been previously created by an extra-heating mecha- 
nism, the central regions of clusters would be characterized by 

an entropy level which does not change with time. Therefore, 

4/3 

after applying the cosmological rescaling, the quantity E z S 
should be positively evolving. However, this evolution turns out 
to be stronger than that measured from the clusters at z > 0.6, 
(1 + z) - 3 , thus excluding that our measured evolution is the ef- 
fect of a z-independent entropy floor. Furthermore, since most 
of our clusters are quite hot systems, they should sample the 
regime where entropy is dominated by gravitational shocks, 
rather than by a pre-collapse entropy floor. Indeed, in the stan- 
dard scenario of pre-heating, one expects a minimum entropy 
level to be created in central regions of groups and poor clus- 
ters, while the ICM entropy structure should be left almost 
unaffected for hot systems, like those included in our sample. 
More in general, the emerging scenario from the entropy prop- 
erties of low- and high-z systems is calling for alternative in- 
terpretations, such as pre-heating with differential entropy am- 
plification, as recently proposed by Voit et al. (2003), or a suit- 
able combination of radiative cooling and heating by feedback 
energy release (e.g., Tornatore et al. 2003; Kay et al. 2003). 

A further possibility is that we are looking at the com- 
bined effect of non-gravitational heating and radiative cool- 
ing. Although each one of these two mechanisms provides an 
excess entropy, hydrodynamical simulations of clusters have 
shown that their combined effect may induce non-trivial con- 
sequences in the entropy properties of the ICM (Tornatore et 
al. 2003). For instance, an episode of impulsive heating may 
cause the creation on an excess entropy, thus increasing the gas 
cooling time and suppressing cooling in central cluster regions. 
Subsequently, cooling starts acting again on the gas, gradually 
decreasing its entropy as gas flows toward the cluster center. In 
this case, one should observe a gradual reduction of the net en- 
tropy detected at the cluster center. Whether this is what we are 
actually observing requires a more sophistical numerical and 
semi-analytical modeling. In any case, there is little doubt that 
the evolution of gas entropy in central cluster regions keeps 
track of the past ICM thermal evolution. As such, its mea- 
surement provides unique information about the latest epoch 
at which a significant amount of energy has been injected into 
the ICM and of the involved energy budget. 

As a diagnostic to connect the thermodynamics of the ICM 
to the past history of star formation, we have investigated the 
relation between entropy and the Iron content. (We remind that 
in these high— z systems Iron is the only element that we can 
sample in their metal budget.) As shown in Fig. ^3 we nave 
a statistically significant indication that values of the central 
entropies are lower in correspondence to metal-richer clusters. 
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Fig. 10. Distribution of the entropy values measured at 0.1 x 
i?2oo versus the Iron abundance measured with respect to the 
solar ratios in Anders & Grevesse (1989). The Spearman's 
rank correlation is significant at 99.0 per cent against the no- 
correlation case. 



This result, along with the evidence that systems at lower tem- 
perature stores lower amount of entropy at the same fraction of 
the virial radius, is an alternative way to show what we reported 
in a previous paper (Tozzi et al. 2003): clusters at lower temper- 
ature appear to have higher Iron abundance (see Baumgartner et 
al. 2003 for a discussion of this trend as observed in the ASCA 
galaxy cluster archive). 

It is now well-established (e.g. Pettini 2003, Finoguenov 
et al. 2002) that most of the star formation and Iron production 
happened at z 2, i.e. before the final collapse of the accreting 
baryons into clusters, and preferentially in low-entropy struc- 
tures like groups along cosmic filaments. Groups are the high- 
z precursor of later-forming rich clusters, which maintain their 
structure and survive for a few crossing times after they merge 
into the cluster environment (as suggested by the frequently ob- 
served cool cores in relaxed systems -e.g. Vikhlinin et al. 2001- 
and corroborated by hydrodynamical simulations -Motl et al. 
2003, and Tormen et al. 2003). This implies that the shorter the 
time elapsed from the aggregation of these clumps in a larger 
system, the larger the number of survival regions at low en- 
tropy and high metallicity. It is worth noticing that these sta- 
ble subclusters should contribute significantly to all the X-ray 
emission-weighted observables because of their high gas den- 
sity. Therefore, from an observational point of view, systems 
with higher metallicity should appear with lower gas tempera- 
ture and be considered as dynamically younger. This supports 
what is shown by the surface brightness profiles, with centrally 
peaked emission, which is characteristic of relaxed systems, 
found to be less evident as more distant clusters are considered. 



We have studied the evolution of the scaling relations between 
observed and derived physical properties of the ICM and dark 
matter potential in a sample of 28 X-ray galaxy clusters with 
redshift in the range z ~ 0.4 — 1.3 (median of 0.73). Using 
Chandra observations, we resolve their surface brightness on 
scales of a few arcsecs (1 arcsec corresponds to between 5 and 
8 /lyQ 1 kpc at these redshifts for the assumed cosmology), out 
to a significant fraction (about 0.7, as median value) of i?soo- 

In general, the surface brightness is well fitted by a sin- 
gle 0— model. A double 0— model is statistically not required, 
suggesting that these distant clusters do not show any signifi- 
cant central brightness excess. This may be indicative that sys- 
tems at high redshift appear as structures in formation and not 
completely relaxed. By collecting all the photons within ~0.5 
xi?500> we have measured a single emission-weighted temper- 
ature (median value of 6.0 keV). The combination of the spatial 
and spectral analysis allows us to recover the gas luminosity 
(observed in the range 1.0 — 116.7 x lO 44 /^ 2 erg s _1 ), gas 
mass (0.6 — 18.1 x 10 13 /i7 5 ^ 2 M Q ) and total gravitating mass 
(1.4 - 21.3 x 10 14 /i 70 1 M Q ) out to i? 5 oo- These quantities are 
used to quantify the evolution of the X-ray scaling relations in 
hot (T gas ^ 3 keV) clusters. 

The slopes of all the investigated correlations tend to be 
higher than what is predicted by self-similar models. The 
largest deviations from the self-similar scaling relations occur 
in the L - T and A/ gas - T relations (Fig. H E EJ. We no- 
tice that the values of the slope propagate through the rela- 
tions in a self-consistent way: for example, the observed de- 
pendence L ~ M 1 ot 88±0 ' 42 and M tot ~ 2 ,1 - 98±0,30 should 
imply L ~ T 3 7±10 that is well in accordance with the mea- 
sured slope of 3.72 ± 0.47. This is in accordance with the ob- 
served properties of gas-dependent quantities, such as luminos- 
ity and gas mass (e.g. Edge & Stewart 1991, Fabian et al. 1994, 
Mushotzky & Scharf 1997, Allen & Fabian 1998, Mohr et al. 
1999, Novicki et al. 2001, Ettori et al. 2002), and suggests that 
simple gravitational collapse is not the only process that gov- 
erns the heating of baryons in the potential well of clusters also 
in the redshift range 0.4—1.3. On the other hand, we do not ob- 
serve any significant dependence of the gas mass fraction upon 
the temperature (see Fig. [HJ, ruling out the simple explanation 
that the increase in the gas budget as function of the mass of the 
system is responsible for the steepening of the L — T relation. 

To explain the observed correlations, a raise in the central 
entropy (and a corresponding suppression of the X-ray emis- 
sivity) is required and can be achieved either by episodes of 
non-gravitational heating due to supernovae and AGN (e.g., 
Evrard & Henry 1991; Cavaliere, Menci & Tozzi 1999; Tozzi 
& Norman 2001; Bialek, Evrard & Mohr 2001; Brighenti & 
Mathews 2001 ; Babul et al. 2002; Borgani et al. 2002), or by se- 
lective removal of low-entropy gas through cooling (e.g. Pearce 
et al. 2001; Voit & Bryan 2001; Wu & Xue 2002), possibly 
regulated by some mechanism supplying energy feedback (e.g. 
the semi-analytical approach proposed by Voit et al. 2002 and 
the numerical simulations discussed by Muanwong et al. 2002, 
Tornatore et al. 2003, Kay et al. 2003). 
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When the slope of the E z M to t — T relation is fixed to 
1.5, its high-z normalization of (0.37 ± 0.02) x 10 u h^M e 
jjgy-3/2 j s we jj m a g reemen t with results from low-z sam- 
ples (e.g. Nevalainen et al. 2000, Allen et al. 2001). This lack 
of evolution follows from the prediction of hydrostatic equi- 
librium within virialized halos. In turn, this is consistent with 
the expectation that the E z M tot -T relation is the least sensi- 
tive to the thermal status of the gas and, therefore, not signifi- 
cantly affected by non-gravitational processes. Still, the M-tot- 
T relation observed at low redshift (e.g., Horner et al. 1999, 
Nevalainen et al. 2000, Finoguenov et al. 2001, Ettori et al. 
2002) is known to have a normalization lower by about 30 
per cent than that obtained from hydrodynamical simulations 
(e.g. Mathiesen & Evrard 2001, cf. also Muanwong et al. 2002, 
Borgani et al. in preparation). It will be interesting to show how 
our measured E z M tot -T for distant clusters compare with pre- 
dictions from hydrodynamical simulations. 

As for the E~ X L - T and E z M gas - T relations, we detect 
a mild negative evolution [B z ~ (—0.1, —1.0), with a devia- 
tion from the non-evolution case B z = by less than 3.2a] 
that becomes more relevant when subsamples with higher cut 
in redshift are considered (see Table0J. For the 16 objects with 
z > 0.6, we measure B z in the range [—1.5, —0.4] with de- 
viation from zero significant at 3.9a. This result points to a 
scenario in which clusters at higher redshift have in general 
lower X-ray luminosities and gas masses than that predicted 
from the self-similar model. This trend appears more clearly 
when the entropy of these hot (T gas <; 3 keV) systems, evalu- 
ated at 0.1i?2oo, is compared with local determinations. A mild 
positive evolution (B z ~ 0.3) is then detected in particular 
when objects at redshift larger than 0.6 are considered. The cut 
at z > 0.6 allows us to improve significantly the fit when a 
X 2 statistic is applied to evaluate the level of evolution. Even 
though a power-law modelization is probably simplistic to rep- 
resent the properties of the ICM as a function of redshift, it 
provides generally a very good fit when clusters with z < 0.6 
are not included. 

As we discuss in Section 4.6, these observed evolutions 
are broadly consistent with predictions from models that pre- 
heat the ICM with a non-evolving entropy floor at the level of 
about 100 keV cm 2 . However, these models are in conflict with 
recent observations of non-isentropic cores in galaxy groups 
(e.g. Ponman et al. 2003, Mushotzky et al. 2003). More de- 
tailed analyses of simulated structure formation are required to 
confirm quantitatively the degree of the expected evolution. 
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RXJ1416: x 2 / dof = 47 - 5 /24 (P<0.01) MS0302: xV dof = 32.8/1 9 (P = 0.03) 
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MS1621: x 2 / dof = 50 - 5 /34 (P = 0.03) RXJ1347: xV dof = 52.4/48 (P = 0.31) 
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Fig. A-l. Data, best-fit /3-model and residuals of the surface brightness profile of the high redshift galaxy clusters in our sample. 
Dotted lines indicate the best-fit background value. The location of R spe c, Ra and Ri a is shown. 
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MS0016: x 2 /dof = 233. 3/167 (P<0. 01) RXJ 1 1 21 : x 2 /dof=39.8/31 (P = 0.13) 
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Fig.A-2. (Continue) 
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MS1137: xVdof= 157.5/1 42 (P = 0.18) 




RXJ1317: x7dof = 64.6/80 (P = 0.90) 
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Fig. A-3. (Continue) 
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RXJ0910: x 2 / dof = 230. 7/193 (P = 0.03) 
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